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Abstract 

The effects of the addition of sub-micrometer sized colloidal silica spheres on 
the linear and nonlinear rheology of semi-dilute solutions of a viscoelastic gel 
are studied. For a 1.4 wt.% solution of the surfactant CTAT, a peak in the zero 
shear rate viscosity rj is observed at approximately equal weight percents of 
silica and CTAT. This peak shifts to lower silica concentrations on increasing 
either the CTAT concentration or the surface charge on silica and disappears 
when the CTAT concentration is increased to 2.6wt%. The increases in ?y 
and the high frequency plateau modulus G on the introduction of SiO"2 are 
explained by considering the increasingly entangled wormlike micelles that are 
formed due to the enhanced screening of the electrostatic interactions. The 
observed decrease in the values of G and r] at higher concentrations of silica 
particles is explained in terms of the formation of surfactant bilayers due to the 
adsorption of the positively charged cetyl trimethylammonium to the negatively 
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charged silica. 

1 Introduction 

The study of the rheology of surfactants, polymers, polyelectrolytes and their 
composites is interesting, not only from an academic point of view, but also in 
the context of possible industrial applications. Viscoelastic surfactant solutions 
are ideal candidates for the study of complex flow phenomena [1] , while mixtures 
of colloids and liquid crystals form novel composite materials under certain con- 
ditions. There has been extensive work on the dramatic changes in the rheology 
of liquid crystals on the addition of colloidal particles [2, 3, 4]. When latex 
microspheres bind to flexible giant vesicles, the resulting membrane distortions 
induce interparticle attractions and the subsequent formation of close-packed 
and ring- like particle aggregates [5]. Recent x-ray studies show that mixtures 
of DNA and cationic or neutral liposomes form stacks of two-dimensional smec- 
tics of aligned DNA chains intercalated by lipid bilaycrs [6] . The experiments of 
Solomon et al. [7] prove that the yielding of concentrated zirconia suspensions is 
extremely sensitive to the headgroup chemical structure of the added surfactant. 
Phase separations in a mixture of hard-sphere rods and colloids have been ob- 
served experimentally, where spontaneous demixing into rod-rich and rod-poor 
states, alternating layers of rods and spheres, columns of spheres in a crystalline 
array etc. are seen to occur [8]. Recent experiments have shown a variety of 
structures (2D hexagonal, lamellar and 2D centered rectangular phase) formed 
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by complexation of cationic and anionic surfactant mixtures with polyclccrolytes 
like DNA, poly-glutonic acid, poly-acrylic acid and polystyrene sulfonate [9]. 

In addition to the systems mentioned above, extensive experimental inves- 
tigations have been cconducted on mixtures of different species of surfactants 
and polymers. Cetyl trimethylammonium bromide (CTAB) micelles show a dra- 
matic increase in length when mixed with a few molecules of oppositely-charged 
gemini surfactants [10]. The results arc explained in terms of the cross-linking 
of the CTAB micelles by the gemini surfactants. Viscometric measurements on 
mixtures of cetyl trimethylammonium tosylate (CTAT) and cetyl trimethylam- 
monium 3-hydroxy naphthalene 2- carboxylate (CTAHNC) show a peak in the 
zero-shear viscosity r/o with increasing HNC~ concentration [11]. The added 
HNC~ is thought to increase the micellar hydrophobicity, which encourages mi- 
cellar growth and increases the solution viscosity. The subsequent decrease in 
r] is understood in terms of the formation of intermicellar connections. Addi- 
tion of the anionic micelle sodium dodecyl sulphate (SDS) to hydrophobically 
modified poly(sodium acrylate) (HMPAA) also shows a peak in the viscosity 
[12]. A maximum in the viscosity of the polymer poly(ethyleneoxide) (PEO) on 
the addition the surfactant SDS is explained in terms of a polymer-surfactant 
complex formation [13]. This peak is shown to correspond to the concentra- 
tion of SDS at which the polymer coils get saturated with the added surfac- 
tant. Addition of mcthyl-Na-benzoates to cetyl pyridinium chloride (CPyCl) 
results in a peak in the zero-shear viscosity [1]. This is explained in terms of 
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a rapid, micellar growth followed by cross-linking of the micelles in the pres- 
ence of the hydrophobic methyl-Na-benzoate. The stress relaxation changes 
from multi-exponential to mono-exponential across the viscosity peak, indicat- 
ing a cross-over from diffusion-controlled to kinetically-controlled processes. It 
is thus easy to control the rheology of polymer and surfactant solutions by 
changing the counterion concentration or by introducing micron-sized colloidal 
additives. Addition of hydrophobic counterions or colloidal particles modifies 
the intcr-micellar interactions and alters the micellar packing considerations, 
and in some cases, may even result in phase transitions and the formation of 
novel complexes. 

In this paper, we discuss our work on the rheology of semi-dilute solutions of 
wormlikc micelles (cetyl trimethylammonium tosylate) in the presence of sub- 
micrometer sized particulate additives (monodisperse silica colloids of diameter 
0.1 /an). In addition to the frequency dependent elastic modulus G'(w) and the 
loss modulus G"(u), we measure the shear rate dependent viscosity 77(7) of the 
samples, and by fitting these to appropriate models, extract the relaxation times 
tr, the zero shear viscosities r) and the high frequency plateau moduli G of 
the mixtures as a function of the weight percent of added silica (Si02) particles. 
We explain our observations by considering the changes in inter-macromolecular 
interactions that arise out of the presence of colloidal additives in the viscoelastic 
gel. 
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2 Experimental 

CTAT powder, purchased from Sigma Chemicals, Bangalore, India, is weighed 
carefully and dissolved at weight percents of 1.4wt.% (31mM), 1.95wt.% (43mM) 
and 2.6wt.% (53mM) in aqueous colloidal silica suspensions prepared at weight 
percents 0.85wt.%, 1.05wt.%, 1.4wt.%, 1.95wt.%, 2.6wt.%, 3.9wt.% and 5.2wt.%. 
The silica colloids were obtained as a gift from Nissan Chemicals, Japan. Pure 
samples of CTAT are also prepared at weight percents 1.4 wt.%, 1.95 wt.% 
and 2.6 wt.% to understand the role of the colloidal silica on the rheology of 
the viscoelastic gel phase of CTAT. In this concentration range, the zero shear 
viscosity rj of CTAT is found to vary with CTAT concentration c as rj <~ c 3 8 . 
The CTAT - silica mixtures are allowed to equilibrate at 30° C for a week. The 
silica suspensions in which CTAT is dissolved are prepared from a stock solution 
of monodisperse colloidal silica (diameter = 0.1 /xm, density = 1.29 g/cc. and 
weight fraction 39.3%), diluted with deionized and distilled water to the required 
concentrations. The effect of electrostatic interactions on micellar aggregation 
in the presence of colloidal Si02 is studied by increasing the surface charge on 
the SiC>2 particles. This is achieved by adding minute quantities of sodium 
hydroxide ([NaOH] ~ 0.178mM) to the silica suspensions. Addition of NaOH 
to Si02 results in the increased dissociation of the surface silanol groups, which 
leads to an increase in the surface charge density of the silica particles [14]. 

All the rhcological measurements are performed in a conc-and-plate geom- 
etry (cone angle = 1°59 , diameter = 4 cm) in an AR-1000N Rheolyst stress- 
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controlled rhcometer (T. A. Instruments, U.K). The sample temperature is con- 
trolled at 25° C for all the experiments. Direct visualization of the samples 
under ambient conditions using optical microscopy shows no evidence of phase 
separation. 

3 Linear rheology 

In viscoclastic gel-forming surfactant solutions, the stress relaxation is governed 
by the breaking and reformation time Ti, rea k and the reptation or curvilinear 
diffusion time r rep of the micelles. In the limit where Tbreak « T rep , the result- 
ing single exponential stress relaxation can be fit to the Maxwell model, where 
the characteristic relaxation time tr = ^/T rep Tb rea k- Fig. 1 shows the frequency 
response data for CTAT 1.4wt.% samples, with (a) Owt.%, (b) 1.3wt.% and (c) 
5.2wt.% silica. The elastic modulus G'(ui) and the viscous modulus G"{lo) do 
not show pure Maxwellian behavior, and are fit to the Cole-Davidson model [15] 
which predicts a stretched exponential form for the stress relaxation and works 
well for many glass- forming liquids. In this model, the dynamic modulus is ex- 
pressed as G*(uj) = G [l — (i + j fa)ra )a ] i where G is the high frequency plateau 
modulus and tr is a measure of an average relaxation time. The inverse of the 
exponent a characterizes the width of the relaxation spectrum and equals 1 for 
Maxwellian relaxation [16, 17]. A decrease in a from its maximum value of 1 
indicates the presence of multiple stress relaxation processes. The frequency 
dependent moduli measured in our experiments are fit to the Cole-Davidson 
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model with Go, tr and a as the fitting parameters. For angular frequencies 
u < 2 rad/sec, the fits (shown as solid lines in fig. 1) are found to improve as 
the concentration of silica additives in the CTAT solutions is increased. The 
deviations observed at high angular frequencies (w > 2 rad/sec) can be at- 
tributed to the effects of breathing (axial stretching modes arising from tube 
length fluctuations) and Rouse modes [18, 19] on the stress relaxation processes 
in the samples. Fig. 2 shows the plot of the fitted values of a as a function of 
increasing SiC>2 content. The values of a, also tabulated in table 1, are found to 
increase monotonically with SiC>2 concentration, implying a separation in the 
relaxation time scales and a tendency towards increasingly mono-exponential 
relaxation on the addition of SiC>2 particles, tr may also be calculated from 
the crossover frequency uu co at which G'(w) and G"(w) are equal, by using the 
relation tr = A good indication of single exponential stress relaxation 

is the semi-circular form of the so-called Cole-Cole figure, where the viscous 
modulus normalized by the high frequency shear modulus G is plotted vs. the 
normalized elastic modulus. Fig. 3 shows the Cole-Cole plot corresponding to 
the data plotted in Fig. 1. The plot for the sample with 5.2 wt.% Si0 2 in CTAT 
approaches the semi-circular form indicative of the dominance of a unique stress 
relaxation mechanism. The Cole-Cole plot for the pure 1.4 wt.% CTAT is almost 
linear and implies a distribution of relaxation time scales, possibly arising out of 
competing processes (for example, comparable breakage and reptation times). 
Figs. 4 (a), (b) and (c) show the plots of the dynamic viscosity rf normalized 
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by the zero-shear viscosity rj (open squares) for the same three samples, where 
rf(uj) = ^ G ' 2 + G "' 2 ._ The dashed lines show the fits to the model for giant worm- 
like micelles given by rf{w) = 2 = [20], where tr is the average relaxation 
time discussed earlier. The fits of rf (ui) to this model are found to work best at 
intermediate SiC>2 concentrations. These results can be interpreted in terms of 
the growth of entangled wormlike micelles, followed by the formation of bilayers 
due to the adsorption of the headgroups on the SiC>2 surfaces. 

4 Nonlinear rheology 

The zero-shear viscosities r] Q of all the samples are also calculated by fitting the 
shear-viscosity (rf(j)) vs. shear rate (7) data obtained in the flow experiments 
to the Giesekus model [21]. The Giesekus model considers the orientation effects 
of flow by introducing a deformation dependent tensorial mobility of breaking 
and reforming micelles and predicts the shear-rate dependent viscosity 77(7) = 
2thV [\/-*- + 4tj?7 2 — 1] . This expression is an analytical solution of a constitutive 
equation for the viscoelastic properties of entangled micelles and assumes an 
adjustable parameter a' — 0.5 [1], where a' is a dimensionless anisotropy factor 
which is related to the relative mobility tensor [3 and the configuration tensor C 
by the relation /? = 1 + ot{C — 1). We would like to note here that experimental 
data for wormlike micellar solutions usually satisfy a' = 0.5 [1]. In addition to 
plots of i]*(ui)/i] vs. oj, fig. 4 also shows the plots of rj(j)/rj (open circles) vs. 
7 for the three samples and the fits to the Giesekus model. A comparison of the 
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values of i] obtained from the fits to the dynamic viscosity and shear viscosity 
data has been made in table 1. The fits to the shear viscosity data are found to 
work best at intermediate SiC>2 concentrations. Fig. 4 may be used to compare 
the degree of validity of the Cox-Merz condition for these samples. The Cox 
Merz rule is an empirical law for wormlike micelles which predicts 77(7) = 
at w = 7. We see that the agreement to the Cox-Merz rule is very poor for the 
pure CTAT sample (fig. 4 (a)) and improves for the intermediate (fig. 4 (b)) 
silica concentration, indicating the growth of giant wormlike micelles. We also 
measure the first normal force difference Ni (7) generated in each sample on the 
imposition of high shear rates and plot the results in Fig. 5. The magnitude 
of the normal stress difference is found to increase with Si02 concentration and 
indicates the formation of structures on the addition of Si02- 

All the experiments and analyses discussed above are repeated for CTAT- 
Si0 2 mixtures with CTAT at concentrations of 1.4 wt.% (+ 0.178 mM of added 
NaOH added to increase the surface charge density of SiC^), 1.95 wt.% and 
2.6 wt.% respectively. The values of tr, G and r) calculated by fitting the 
data to the models discussed above are plotted in Figs. 6-10. In fig. 6, we 
plot the results for CTAT 1.4wt.% + Si02- In fig. 7, we compare the values of 
t r , Go and r/ for CTAT 1.4wt.% + Si0 2 , without added NaOH (solid circles) 
and with 0.178mM NaOH (open triangles). All the plots for the mixtures with 
untreated Si02 show prominent peaks at Si02 weight percents between 1.3 and 
1.9. The peaks clearly shift to lower Si02 concentrations when Si02 is treated 
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Table 1: The parameters a, G , tr and rj , obtained by fitting the CTAT 1.4 
wt.% + SiC>2 data to the models described above, a and G are obtained from 
fits to the Cole Davidson model, the mean relaxation time tr is the inverse 
of the crossover frequency loco" 1 and rj is obtained by fitting to the dynamic 
viscosity (DV) and the Giesekus (GM) models respectively. 
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with NaOH. Figs. 8 and 9 are the plots of tr, G and r] for CTAT-silica 
mixtures with the CTAT weight percent fixed at 1.95 and 2.6 respectively. In 
table 1, we have compiled the values of the fitting parameters obtained from 
the fits to the frequency response and steady shear experiments on CTAT 1.4 
wt.% as a function of added Si02 concentration. The values of r) obtained from 
fits to the dynamic viscosity and the Giesekus models are found to agree more 
at intermediate Si02 concentrations, indicating a better agreement with the 
empirical Cox-Merz rule under these conditions. In fig. 10, we have plotted r/ , 
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and tr for the different CTAT concentrations. For CTAT 1.4wt.%, the values 
of r]o, Go (not shown) and tr arc found to peak at silica concentrations ~ 
CTAT concentrations. The peaks shift to lower SiC-2 concentrations for CTAT 
1.95wt.%. For CTAT 2.6wt.%, the peaks disappear or possibly shift to very low 
silica concentrations that lie beyond the range of the present measurements. 

5 Discussions 

As discussed earlier, the stress relaxation in viscoelastic solutions of worm- 
like micelles is dominated by two processes, reptation and reversible scission 
[1, 16, 17]. The value of a in the Cole-Davidson model [15] gives an estimate of 
the width of the distribution of relaxation times. Multi-exponential relaxation 
processes are seen in dilute micellar solutions where a < 1, with a increasing 
to 1 with increasing surfactant and/or salt concentrations [22, 23]. On adding 
SiC-2 to semi-dilute aqueous solutions of CTAT, we find a similar increase in the 
values of a suggesting an approach towards single exponential stress relaxation. 
From our viscosity data, it is clear that there is an increased tendency towards 
behavior characteristic of giant wormlike micelles, followed by significant depar- 
tures on further addition of Si02 particles. 

In this section, we explain the changes in the structure and dynamics of 
the cylindrical CTAT micelles when silica colloids are added to CTAT solutions 
under appropriate conditions. In the case of CTAT 1.4 wt.%, an addition of 1.3 
wt.% Si02 particles results in the relaxation time increasing by 600%, G Q in- 
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creasing by almost 37% and ri increasing by 1600% as compared to pure CTAT 
solutions. We find that the peaks in r] and tr shift to lower concentrations of 
SiC>2 as the surfactant concentration is increased from 1.4 to 1.95 wt.%, and dis- 
appear for CTAT 2.6wt.% (fig. 10). The increased fluidity of the samples on the 
addition of SiC-2 can be explained by considering the adsorption of the surfac- 
tant headgroups to the silica surfaces due to the attractive interactions between 
them [24]. The surfactants can be thought to form bilayers around the silica 
surface, with the headgroups constituting the outer layer and with the silica at 
the center. Enhanced formation of bilayers with increasing SiC>2 concentration 
can therefore explain the observed decreases in 770, Go and tr. At lower SiC>2 
concentrations, giant wormlike micelles outnumber the bilayers formed, and the 
presence of Si0 2 in solution serves to screen the electrostatic interaction between 
micelles, resulting in increased micellar entanglement, which manifests as peaks 
in plots of the relaxation times, shear moduli and viscosities of the samples. 
Addition of NaOH to 1.4wt.% CTAT results in a decrease in the magnitudes 
of rj , G and tr, with the peak shifting to lower Si02 concentrations (fig. 7). 
Treating the Si0 2 particles with NaOH increases their surface charge density by 
increasing the dissociation of the surface silanol groups. Increasing the surface 
charge on silica causes them to bind more strongly to the CTA + , which results 
in the formation of bilayers at lower silica concentrations, clearly indicated by a 
shift in the peaks to the left (fig. 7). The dramatic changes in the intermicellar 
interactions on treating the SiC-2 particles with NaOH are thus responsible for 
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the observed changes in tr, Go and r/ . On continued addition of silica to the 
CTAT sample, a stage comes when all the headgroups get adsorbed onto the 
silica surfaces, leaving very few entangled micelles in solution. 

Next, we focus on the shift in the viscosity peak to lower concentrations of 
added silica, when the surfactant concentration is increased to 1.95wt.% (fig. 8) 
and 2.6 wt.% (fig. 9). It is useful to mention here that these weight percents lie in 
the semi-dilute concentration regime of CTAT [25] . An increase in the surfactant 
concentration is known to result in an increase in the degree of ionization of 
the micelles [19]. Higher surfactant concentrations encourage the formation of 
bilayers at lower silica concentrations, resulting in the shifts of the peaks to lower 
silica concentrations. We conclude this section by stressing that we identify the 
peak as a signature of the formation of longer "worms" on the addition of SiC-2. 
We believe that this is followed by the formation of bilayers, resulting from the 
subsequent adsorption of the micellar headgroups to silica, which is manifested 
experimentally as a decrease in tr, G and i] . 

6 Conclusions 

In this paper, we study the modifications in the rheology of semi-dilute solu- 
tions of CTAT on the addition of submicrometer sized silica spheres. The drastic 
changes in rheology that we observe as a function of increasing silica concen- 
tration are explained by considering the electrostatic interactions between the 
surfactant and silica particles. We would like to note here that modifications 
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of the rheology of surfactant solutions by the controlled addition of particulate 
matter can have diverse uses in the industry. 
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Figure 1: The clastic modulus G'(ui) and the viscous modulus G"(u>) vs. an- 
gular frequency lo for CTAT 1.4wt.%+Si0 2 samples. The solid lines show the 
fits to the real and imaginary parts of the Cole-Davidson model. The SiC>2 
concentrations corresponding to each frequency response curve are (a) wt.%, 

(b) 1.3wt.% and (c) 5.2wt.%. 

Figure 2: Plot of the values of a vs. S1O2 concentration for CTAT 1.4wt.% 
solutions, a is obtained by fitting the frequency response curves to the Cole- 
Davidson model. 

Figure 3: Cole-Cole plots for CTAT 1.4 wt.% with Si0 2 concentrations wt.% 
(open squares), 1.3 wt.% (plus-centered circles) and 5.2 wt.% (filled diamonds) 
respectively. 

Figure 4: The normalized dynamic viscosity i]*(uj)/rj vs. lo (open squares) and 
the normalized shear viscosity r)(j)/rj vs. shear rate (open circles) 7 for CTAT 
1.4wt.%+Si02 samples. The dashed lines show the fits of the dynamic viscosity 
to ii* = , n ° „ and the solid lines show the fits to the Giesekus model. The 

Si02 concentration corresponding to each plot is (a) wt.%, (b) 1.3wt.% and 

(c) 5.2wt.%. 

Figure 5: The first normal stress difference N\ vs. shear rate 7 for CTAT 
1.4wt.%+Si0 2 samples. 

Figure 6: Plot of viscoelastic parameters for CTAT 1.4wt.%. (a) shows the 
calculation of the average relaxation time tr from the crossover frequency (cross- 
centred triangles), fits to the dynamic viscosity model (open circles) and the fits 
to the Giesekus model (solid circles) respectively, (b) shows the plots of the 
values of G at different Si02 concentrations obtained from the fits to the Cole- 
Davidson model, (c) shows the values of r/ obtained from the fits to the dynamic 
viscosity model (open diamonds) and the Giesekus model (solid circles). 

Figure 7: Plots of the fitted parameters in the absence (filled circles) and in the 
presence (open triangles) of 0.178mM NaOH for CTAT 1.4wt.%. (a) and (b) 
show the calculation of the zero-shear viscosity r) from the fits to the Giesekus 
and the dynamic viscosity models respectively and (c) shows the plots of the 
values of G Q at different Si02 concentrations obtained from the Cole-Davidson 
fits, (d) shows the values of tr obtained by using the relation tr <~ lo^ . 
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Figure 8: Plots of the fitted parameters for CTAT 1.95wt.% + Si0 2 . (a) shows 
the plot of the average relaxation time tr obtained from the inverse of the 
crossover frequency (solid triangles), fits to the dynamic viscosity model (open 
diamonds) and the fits to the Gicsekus model (solid circles), (b) shows the 
plots of the values of G Q at different SiC>2 concentrations obtained from the 
Cole-Davidson fits and (c) shows the values of r] obtained from the fits to the 
functional form suggested for 77* (open diamonds) and the Giesekus model (filled 
circles). 



Figure 9: Plots of the fitted parameters for CTAT 2.6wt.% + Si0 2 . (a) shows the 
plot of the average relaxation time tr i) obtained using the relation tr = uj^ 1 
(solid triangles), ii) from fits to the dynamic viscosity model (open diamonds) 
and (iii) from the fits to the Giesekus model (solid circles), (b) shows a plot of the 
values of G at different Si02 concentrations, obtained from the Cole-Davidson 
fits and (c) shows the values of r) obtained from the fits to the dynamic viscosity 
model (open diamonds) and the Giesekus model (solid circles). 



Figure 10: Plot of the zero-shear viscosity r] obtained by fitting to the dynamic 
viscosity data for of 1.4wt.% (circles), 1.95wt.% (up-triangles) and 2.6wt.% 
(down-trianglcs) CTAT. The relaxation times tr obtained from u co are also 
plotted vs. silica concentration for the three concentrations of CTAT and are 
indicated by the same symbols described above. The dotted lines are guides to 
the eye. 
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This figure "rbfig2.PNG" is available in "PNG" format from: 



http://arXiv.org/ps/cond-mat/0409435vl 
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